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Abstract 

We study density perturbations, including their non-Gaussianity, in models in 
which the decay rate of the curvaton depends on another light scalar field, denoted 
the modulaton. Although this model shares some similarities with the standard cur- 
vaton and modulated reheating scenarios, it exhibits interesting predictions for /nl 
and ^NL that are specific to this model. We also discuss the possibility that both 
modulaton and curvaton fluctuations contribute to the final curvature perturbation. 
Our results naturally include the standard curvaton and modulated reheating sce- 
narios as specific limits and are thus useful to present a unified treatment of these 
models and their variants. 



1 Introduction 



Cosmological observations are increasingly precise, providing us with a lot of information 
on the origin of cosmic structure, i.e. on primordial density fluctuations. Although the 
fluctuations of the inflaton field are considered as the main candidate for their origin, 
other possibilities have also been investigated (see e.g. [1] for introductory lectures), es- 
pecially in the light of recent constraints on primordial non-Gaussianity. The degree of 
non-Gaussianity in primordial density fluctuations is usually characterized by the non- 
linearity parameter /nl which represents the amplitude of the bispectrum. In the case 
of standard slow-roll single-field inflation, /nl is predicted to be too small to be observ- 
able. On the other hand, the present constraints on /nl for local type non-Gaussianity 
obtained from cosmic microwave background (CMB) and large scale structure are respec- 
tively fSl""-^^ = 37.2 ± 19.9 (la C.L.) from WMAP9 [2] and /jj""'"'^ = 62 ± 27 (la C.L.) 
from NRAO VL A Sky Survey [3] , which may give some hints that the value of /nl is away 
from zero. 

In this context, other candidates for primordial fluctuations, in particular those giving 
significant /nl, have been extensively discussed such as the curvaton model [4-6], mod- 
ulated reheating scenario [7,8], inhomogeneous end of hybrid inflation [9-12], modulated 
trapping [13, 14] and so on. Even if we limit ourselves to the curvaton and modulated re- 
heating scenarios, various extensions of them have been proposed and studied, for example 
mixed inflaton-curvaton model [15-19], mixed inflaton- modulated reheating model [20], 
multi-curvaton [21,22], modulated curvaton [23,24] and so on*^. In most of those scenar- 
ios, a light scalar field (degree of freedom) other than the inflaton is involved in some way, 
and the final values of density perturbations depend on how the initial fluctuations are 
converted to the final ones during the evolution of the early Universe. This consideration 
provides a strong motivation to treat models involving curvatons and/or modulatons in 
a unified formalism (see [34]). In the present work, we wish to focus our attention on a 
scenario that interpolates between the modulated reheating and the curvaton models: that 
of the modulated curvaton decay, in which the decay of the curvaton field is modulated 
by the dependence of its decay rate T on another fluctuating scalar field, which we will 
call the modulaton*^. 

In this paper, we focus on this new mechanism of generating primordial density fluctu- 
ations and derive its predictions, paying particular attention to non-Gaussianities by com- 
puting the bispectrum and trispectrum. Although the model we propose here is, in some 
sense, a straightforward extension of the curvaton and modulated reheating scenarios, we 
find that it leads to a rich phenomenology, with interesting observational implications for 

"^^In most works, the curvaton potential is assumed to have a quadratic form, however, other types of 
the potential have also been discussed such as curvaton with self-interaction [25-29] and pseudo-Nambu- 
Goldstone curvaton [30-33]. 

#^Our scenario differs from the so-called modulated curvaton model, discussed in [23,24], where the 
curvaton first plays the role of a modulaton during the decay of the inflaton, then decays at some later 
time as in the usual curvaton model. 
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primordial non-Gaussianity. We also provide general formulas, which could be applied to 
other similar types of scenarios (see [34] for a systematic approach, including isocurvature 
perturbations) . 

The structure of this paper is as follows. In the next section, we derive a general 
expression describing the final curvature perturbation, up to third order in the perturba- 
tions. Focussing then, in Section 3, on the modulaton fluctuations, we analyse the density 
perturbation and compute the non-linear parameters such as /nl and (^nl- In Section 4, we 
investigate models in which the curvaton fluctuations also contribute to the observed per- 
turbations in addition to those from the modulaton. The final section contains a summary 
of this paper. 

2 Computing the post-decay perturbation 

The present scenario relies on the presence of three scalar fields: an infiaton field (f) which 
drives the inflationary expansion; a curvaton (or modulus) field a, with an energy density 
negligible during inflation, which, well after inflation, oscillates at the bottom of its po- 
tential before decaying; and, finally, a modulaton field Xi which is light during inflation 
and thus acquires fluctuations from the amplification of quantum fluctuations. The crucial 
assumption here is that the decay rate of the curvaton a depends on the modulaton x- 
Therefore, fluctuations of x directly lead to a varying decay rate and eventually produce 
density fluctuations. 

In practice, we will not need any detail about the infiaton field. Its role will be simply to 
drive inflation so that the modulaton field can acquire some fluctuations. In our scenario, 
a can be either light during inflation (mo- <^ H), in which case it will also acquire some 
fluctuations, or be massive (mg- ^ H) in which case its fluctuations are suppressed. 
Strictly speaking, the curvaton scenario assumes a light field during inflation, but there 
exist models where the additional scalar fields, such as moduli, are not necessary light 
during inflation. Whereas the original curvaton scenario would not apply to scalar fields 
of this type, our model does. In the following, although a is a scalar field, we will be 
interested in the cosmological phase where it oscillates at the bottom of its potential and 
can be effectively described as a fluid, which is pressureless if the potential is quadratic. 
Note that our formalism also applies to the decay of the infiaton oscillating in a quadratic 
potential at the end of inflation, if a is simply replaced by the infiaton 0. Our formalism 
thus includes automatically the modulated reheating scenario. 

For each fluid characterized by an equation of state parameter Wi = Pi/pi, which 
is assumed here to be constant, it is convenient to introduce the non-linear curvature 
perturbation Q [35] (see also [36-38] for a covariant definition) 



where 6N denotes the local perturbation of the number of e-folds and a barred quantity 
must be understood as homogeneous. From the above formula, the nonlinear energy 
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density of the species i can be written locally as 

p,(t,x) = p,(t)e3(^+"'>)(^-^^). (2) 

In our case, we will consider only two species: radiation {wr = 1/3) and the curvaton field, 
treated pressureless fiuid = 0). 

Using the instantaneous decay approximation, the value of the Hubble parameter at 
the decay of the curvaton a (or, alternatively, of the infiaton to describe inhomogeneous 
reheating) is given by 

H = T{x): (3) 

where the decay rate F is a function of the modulaton x- Because of the modulaton 
fiuctuations 6x = H/ (27r) generated during infiation, the decay hypersurface characterized 
by the above relation is inhomogeneous. Using Friedmann's equations*'^, this implies for 
the local energy density 

PtotitD) ~ H^tn) ~ f2(t^) ' ^ ^ 

where t^ix) represents the local decay time. Substituting 

PtotitD, x) = J2p^itD,x) = ^p,(tz5)e3(^+'"')K-^^-) (5) 

i i 

in the relation (4), both for the matter contents just before and just after decay, we find 

(1 - fi.) e^(^'-'^^-) + e^K-^^-) = (1 + 5vf = e^^"^^-), (6) 

where we have introduced the curvaton fraction of the total energy density (just before 
the decay) VL„ = Pajptot, as well as the (nonlinear) relative fiuctuations of the decay rate 

_ V{tD,x) 

The first equality in (6) gives us the expression of 6N£, as a function of the two pre-decay 
curvature perturbations Cr and Co-- And the second equality in (6) yields the expression 
of the post-decay curvature perturbation ( (carried by the only-remaining radiation fiuid) 
function of SNd, namely 

( = 6ND + ^Hl + 6r). (8) 

There is no general nonlinear expression for 6N£) given in terms of Cr and Co-, but by 
expanding the first equality of (6) order by order in the perturbations, one can iteratively 



"^^Note that we are implicitly using the separate Universe approach "where distinct regions are described 
by FLRW universes. This is justified by the fact that the perturbations we are interested in are on 
super-Hubble scales at the time of the decay. 
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obtain an explicit expression for ^A^^^ valid up to any order. Computing SNn up to third 
order with this method and substituting in (8), we finally get the following expression for 
the post-decay curvature perturbation: 

C = Cr - ^Sr + ^5 - 4 \ir^ + 2r-9)6l + 4: (r^ + 2r - 3) 5' - 4 (r^ + 2r - 3) S6r] 
b o 72 

+ [(-3r^ - lOr^ + 22r^ + 54r - 135) 6^ + 8 (3r^ + lOr^ - Ar^ - 18r + 9) 
1296 

-12 (3r^ + lOr^ - Ar^ - 18r + 9) S^6r + 6 (3r^ + lOr^ - lOr^ - 30r + 27) S6l] , 

(9) 

where we have introduced, for convenience, the curvaton isocurvature perturbation 

S = 3(C. - Cr) (10) 

and the parameter r, defined by 



Apr + 3p^ 



Although isocurvature fluctuations can also be generated in principle, we restrict our 
analysis to adiabatic perturbations in the present work (see [34] for an analysis including 
isocurvature modes). 

Note that the perturbations Cr, Co- and 5y are related to the fluctuations of the inflaton, 
curvaton and modulaton, via the expressions'^^ 

(p a 3 a-^ F 2r 6r 

where a prime denotes the derivatives with respect to x- Substituting the above expressions 
in (9) would thus give the curvature perturbation C as a function of the fluctuations S(f), 
Sa and Sx- 

Once the curvature perturbation has been computed, here up to third order, it is useful, 
in order to confront the model with observations, to calculate the power spectrum P^, 
bispectrum and trispectrum T^^. They correspond, respectively, to the 2-point, 3-point 
and 4-point correlation functions in Fourier space: 

(CfciCfc.) = {27rf 6{k, + k2)Pdki) , (13) 
{Ck.CUks) = (27r)35(fci + k2 + h)B^{h, /C2, h) , (14) 
{CkMk-Xk,) = (27r)3(5(fci + ^2 + ^3 + ki)T^{h, k2, h, h) . (15) 



"^^We include only the linear term for the inflaton fluctuations, since their non-Gaussianity can be 
neglected in the simplest models. 
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In the case of local non-Gaussianity, it is convenient to express the bispectrum and trispec- 
trum in terms of the power spectrum and to introduce the so-called non-linearity param- 
eters /nl for the bispectrum, tnl and (7nl for the trispectrum: 

Bcih, h, h) = ^/nl {Pdh)Pc{h) + P^{k2)Pdh) + p^{h)P^{h)) , (16) 

T^(/ci, k2, h, ki) = Tnl (Pc(^i3)-Pc(^3)-Pc(^4) + H perms.) 
54 

+^9nl {Pc;{k2)Pdh)Pdh) + 3 perms.) . (17) 
Quite generically, if the curvature perturbation can be written in the form 

C = NJif- + ^Nab Sif^'Sif' + ]:Kbc Sif^'Sif'Sif'^ + ■■■ , (18) 
2 

where the ip'^ denotes any number of light scalar fields, labelled by the index a, with 
statistical independent fluctuations generated during inflation*^, 

mM.) = i^^)'^ik, + k2)P.{k,)6^'', P^{k,) = ^V., T^*={§^) ' (19) 
the power spectrum is given by 

P^ = NaN''P,, (20) 
and the non-linearity parameters by the simple expressions [39-41] 



6 ^ NaNhN-' 
5 {N^N" 



/nl = -^/nl — -7T7T7~^^ (21) 



TNL = — 22 

(iVdiV^)' 
. _ 54 iV,fe,iV«iV''iV'= 
= 25^^^ - (iV,iV.)3 ' 

where we have used the Kronecker symbols to raise the scalar field indices, e.g. iV" = 6°'^Nb. 

Equipped with the general formalism presented in this section, we now discuss in more 
details various scenarios based on the modulated decay of the curvaton in the following 
sections. 



3 Modulaton dominated case 

In this section, we restrict ourselves to the simplest case where only the modulaton fluc- 
tuations are relevant, i.e. we neglect the fluctuations of the curvaton and of the pre-decay 

"^^Wc implicitly assume, for simplicity, slow-roll inflation with standard kinetic terms for all scalar fields, 
so that their fluctuations are approximately Gaussian. 



5 



radiation fluid. Note that the curvaton fluctuations can be ignored in two distinct situa- 
tions: firstly, if their contribution turns out to be numerically negligible in the expression 
for (] secondly, if the curvaton field is massive during inflation (i.e. rricr ^ H), in which 
case its fluctuations are suppressed. 

From the general formula (9), the curvature perturbation in this case reduces to 

r (3r^ + lOr^ - 22r'^ - 54r + 135) 



C 



6 



72 



(r^ + 2r - 9) 5^ - 



1296 



(24) 



Substituting the expansion of 6r in terms of 6x, given in (12), one obtains, up to third 
order, the expression 



c 



r r' r 
"gY ^~ 72 



6- + {r' + 2r-9) — 



1296 



+ (3r^ + lOr^ - 22r^ - 54r + 135) 



5x' 



r'" rT" 

36— + 18 (r^ + 2r - 9) -— 
1 1 ^ 

(25) 



As mentioned earlier, our analysis includes the original modulated reheating scenario if a 
is identified with the inflaton. One can indeed check that the above formula, in the limit 
r = 1 (since the inflaton dominates the total energy density), reduces to the third-order 
expression obtained for ( in the modulated reheating scenario [42]. 

Since the expression (25) is exactly of the form (18) with x as a unique scalar field, 
one can readily use the general formulas (20-22) to determine the power spectrum and the 
non-linearity parameters. The power spectrum is thus 



r \ 2 

6 



^1 P. 



while we obtain for the non-linearity parameters 



NL 



9nl 



- 3-2- 



Y'2 



- 2 



T"NL = /] 



2 

NL ; 



(26) 



(27) 



36- 



r^r' 



r'3 



18 (r^ + 2r - 9) 



rr" 

p/2 



3^4 ^ ^Q^d _ 22r' - 54r + 135 



(2^ 



In the limit r = 1, one recovers the predictions of the modulated reheating scenario. In 
the limit r <^ 1, one finds that the dominant terms in /nl and g'NL; barring some special 
cancellation, scale like 1/r and 1/r^, respectively. This is different from the standard 
curvaton model, where both /nl and (^nl scale like 1/r*^. This specific feature of our 
model can be expressed by the relation 



5'NL 



^/nL' 



(r < 1) 



(29) 



T^^This is true only if the curvaton potential is quadratic, which we assume in the present work. When 
the potential is not quadratic, gNL oc l/r^ [26]. 
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where the numerical value of the coefficient C, and in particular its sign, depends on the 
functional form of F: 

C^I{^,./-^-^sI^){-2nl,sy\ (30, 

To go further and determine quantitatively the parameters /nl and (7nl, taking into 
account the observed amplitude of the power spectrum, we need to assume specific ex- 
pressions for the function T{x)- We consider below three possibilities for r(x). 



Case I:r(x) = roxMP>0). 

By imposing the CMB normalization for the power spectrum V( = 2.4 x 10~^ [2], 
Eq. (26) determines the parameter r in terms of the inflationary Hubble parameter Hi^f 
and X- 

1-8 X 10^3 X .... 

Furthermore, r should be less than unity by definition, which limits the parameter range 
of -ffinf and X to satisfy 

4.5 X 10^^ GeV x .... 

Regarding the non-linearity parameters /nl and s'nl, the dependence on the functional 
form of r only appears in the combinations rr /(r')^andrr /(F')^. By substituting the 
functional form, we obtain 

/,, = !(|^-2-r. (33) 

pr 



„ _ p2 (3^4 + I0r3 - 4r2 - 18r + 9) - 18p (r^ + 2r - 3) + 72 
5'NL — ■ 



(34) 



In the limit r <C 1, one thus finds that /nl is always positive and can become large if r 
is small enough. In the same limit, (7nl is positive and enhanced by the factor as 
already mentioned. In Fig. 1, we show contours of /nl (left plot) and gm. (right plot) in 
the x^-f^inf plane. The value of r is fixed by the CMB normalization as described above. 



Case II: F(x) = Fq 



M \mJ 



In many models, the coupling can be written as a Taylor expansion of x? where M 
represents some high energy scale and x/^ 1 is assumed. The coefficients a and b are 
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Figure 1: Contours of /nl (left) and (7nl (right) in the x^-^^mf plane for the Case I. p = 2 
is assumed. Shaded regions correspond to the parameter space where we cannot obtain 
the right amplitude for P^. 



parameters of order one, which are supposed to depend on the details of some explicit 
model of high energy physics. The non-linearity parameters are then given by 



126 



NL 



fi'NL 



?-2. 

r 



366 



(r^ + 2r - 9) + 3r^ + lOr^ - 22r^ - 54r + 135 



(35) 
(36) 



where we have used ^ 1- As easily read off from the above expressions, the signs of 
/nl and (/nl can be positive or negative, depending on a and b. Since a and b are assumed 
to be C(l), the amplitude of /nl and (/nl is mainly controlled by the value of r. 



Case III: r(x) = To 



M 



{q > 2). 



In the latter case, we consider the possibility that the expansion starts with a higher 
order polynomial. Once again, M is some high energy scale characterizing the underlying 
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physics, and we assume xl^ 1- The non-hnearity parameters are now given by 



NL 



6 
5 

54 1 



6g- 1 (M\^ 



q \X J 



9 

r + --2, 

r 



36 



(g-l)(g-2) (M\ 



2q 



- 54r 



+ 18 {r 
135)] . 



(37) 



2r - 9) 



(3^ 



Since xl^ ^ 1; oiis sees that /nl can be large even with r = 1, but the sign of /nl is 
then negative, which is in contradiction with the current constraints from WMAP [2]. On 
the other hand, when r -C 1, /nl and (^nl can be both positive, which is similar to the 
case I. 



4 Hybrid curvaton-modulaton case 

In this section, we consider the general situation where both the curvaton and modulaton 
fluctuations contribute to the final density perturbation. We also take into account the 
inflaton contribution in the power spectrum. 

The curvature perturbation in this case is given by (9), with the substitutions of (12). 
This leads to an expression of the form (18), which contains now three scalar fields </>, a 
and X- If is convenient to introduce two dimensionless parameters that characterize the 
relative contributions of a and x fo the total power spectrum, defined by 

iV2 Nl 
- - " - ^ (39) 



The inflaton contribution in the power spectrum is thus given by 1 — Ho- — H^. 

The bispectrum parameter /nl can then be decomposed into three terms (the inflaton 
does not contribute to the non-Gaussianity here) 

/nl = -a /nl ^ + 2 "^a '^x f^L^ + ^NL ^ (^0) 



where we have defined 



/NL ^2' /NL /NL ^2 

Their explicit expressions are 



(42) 

= -r + l-2, (43) 
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where the first equation corresponds to the standard curvaton expression [39], while the 
second one is the modulaton contribution calculated in the previous section. The final one 
is a mixed contribution. 

Similarly, the trispectrum coefficient ^nl can be decomposed into four terms: 



where 



~(<73) 

9nl 



9nl 



9nl 



N. 



~(x^) 



iV3 ' 

By using (9) and (12), we can explicitly write down these quantities as 
3r^ 



ATS 

X 



9nl 



9 1 
lOr - - + 



~9^.^y 



9^L ^ 



3r' 
1 



9 15 5 
+ lOr -, 



2r 



2r - 3) ((3r2 + 4r - 9) + ^ (r^ + 2r - 3) 



p/2 ' 

36^^ + 18 (r^ + 2r - 9) + (3r^ + lOr^ - 22r2 - 54r + 135) 



p/3 



(44) 
(45) 

(46) 
(47) 
(48) 
,(49) 



where one can identify the usual curvaton contribution [43], the pure modulaton con- 
tribution calculated in the previous section, as well as two mixed curvaton-modulaton 
contributions. 

Finally, the tnl coefficient can be decomposed as 



tnl 



'NL 



+ 2 ^o- ^X /nL Vnl''' + ^cr ('^o- + '^x) ( /] 

I 9^ ^2 fix^) I ^3 f(x^)^^ 

-t-^ ^o- --^ JnL iNL "T "X I •'NL J ■ 



c(o-x) 
'NL 



(50) 



In the r <^ 1 limit, one finds that the dominant terms for the non-Gaussianity coeffi- 
cients are 



NL 



1 

2r 



3El + 12S.S^ 



18S^ - 12S^ 



rr" 

X p/2 



(51) 



Tnl 



4^2 



^S^S^ + (28-16 



rr 

r2 



Et 36 - 48- 



rr" 



2r 2 



+ 16 



r^r 



S'NL 



9 - 

^ ^X 

2r2 



3 (6S.S, + + lOS^) - 12S^ (S. + 3S, 



rr" 



+ 8S 



p'4 
■p2-p///" 



oil 



r 
(52) 



"X/ -1-^/2 ' """X p/3 

^ '(53) 

Here it is interesting to notice that the leading term of (^nl vanishes when = 0. Thus 
the enhancement by the factor of l/r^ comes from the modulaton fiuctuations, which is 
absent in the standard curvaton model. This property of the trispectrum could be useful 
to discriminate this model from other ones. 
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5 Summary 



We have investigated the density perturbations in a scenario based on the modulated decay 
of the curvaton, where the curvature perturbation is generated via the modulation of the 
curvaton decay rate due to its dependence on another scalar field, called the modulaton x- 
We have paid special attention to non-Gaussianity, since it is potentially the best way to 
discriminate among various scenarios, and we have computed specifically the non-linearity 
parameters /nl and (^nl in this class of models. 

As discussed in Section 3, this model shares some similarities with the usual curvaton 
and modulated reheating models: the size of these parameters are mainly determined by 
r, which is similar to the curvaton model. On the other hand, the signs of these parame- 
ters highly depend on the functional form of r(x), which is the same as usual modulated 
reheating model. However, the model also exhibits interesting predictions coming from 
a hybrid nature of this model. When the curvaton is a subdominant component of the 
Universe at its decay, the non-linearity parameters are related as (^nl oc /^l where the 
proportionality factor depends on the functional form of r(x). Thus this model generi- 
cally predicts enhanced (^nl? which is different from the standard curvaton and modulated 
reheating models, predicting g^i^ ~ /nl- 

We have also investigated a more general case where fluctuations of the curvaton itself 
also contribute to density fluctuations. We have presented the formulas of the curvature 
perturbation up to the 3rd order, and the non-linearity parameters. The formulas given in 
such a case naturally also include the standard curvaton and modulated reheating models, 
which provides a unified treatment of those kind of models and their variants. 

Although we have considered only adiabatic perturbations in the present work, the 
modulated decay of the curvaton could also produce isocurvature perturbations, similarly 
to the standard curvaton scenario [44]. It would be interesting to study these possible 
isocurvature modes and their non-Gaussianities following the analysis introduced in [45, 
46]. Isocurvature modes, possibly correlated with adiabatic modes, lead to very specific 
signatures in the CMB non-Gaussianities and Planck or future CMB data could detect 
these isocurvature non-Gaussianities [47,48]*'^. 

In the near future, one can hope that new cosmological data of unprecedented precision, 
in particular from Planck, will enable to test the model presented here, together with 
various other mechanisms that generate primordial density perturbations. In this respect, 
the unified treatment proposed in this paper should be useful for a simplified confrontation 
of a large class of models with cosmological data. 

Note added: While completing this manuscript, we became aware that the authors 
of [52] were working on a very similar topic. 

'^''For constraints on isocurvature non-Gaussianities from current data, see [49-51]. 
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